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Abstract 

We present explicit solutions for the ordinary differential equations system describing 
the motion of the particles beneath small-amplitude capillary-gravity waves which 
propagate on the surface of an irrotational water flow with a flat bottom. The required 
computations involve elliptic integrals of first kind, the Legendre normal form and a 
solvable Abel differential equation of the second kind. 

1 Introduction 

We consider the problem of water waves in a domain of finite depth bounded 
above by a free surface and under the combined effects of gravity and surface 
tension. We suppose that the water flow is irrotational. Mathematically, the 
problem is formulated as a free boundary problem for incompressible Euler 
equations with the irrotational condition. After rewriting the equations in an 
appropriate non-dimensional form, we have two non-dimensional parameters S 
and e, the shallowness parameter and the amplitude parameter, respectively, 
and another non-dimensional parameter W e called Weber number, which comes 
from the surface tension on the free surface. We simplify the governing equa- 
tions with a linearization which is slightly different from the classical case in 
line with the Stokes condition for irrotational flows (see, for example, [5], [TU]V 
By this linearization, we obtained a parameter cq by which we can describe 
different backward flows in the irrotational case: still water (c = 0), favorable 
uniform current c > 0, adverse uniform current Cq < 0. 

Further, we get the general solution of the linearized problem. Notice that there 
are only a few explicit solutions to the nonlinear governing equations: for grav- 
ity water waves, Gerstner's solutior£j][l6] and the edge wave solution related to it 



1 This solution was independently re-discovered later by Rankine 1301 . Modern detailed 
descriptions of this wave are given in the recent papers [2] and |20| . 
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(see [3]), for capillary water waves, Crapper's solution [T3] and its generalization 
in the case of finite depth (see [26]). 

After getting the general solution of the linearized problem we investigate the 
nonlinear equations of the motion of the fluid particles. In the case the con- 
stant Co equals the non-dimensional speed of propagation of the linear wave, 
the required computations involve elliptic integrals of first kind and their Leg- 
endre's normal form. The exact solutions obtained in this case contain in their 
expressions Jacobian elliptic functions. Only one solution is presented in detail, 
the others will be presented in a future paper. In the case the constant Co is 
different from the non-dimensional speed of propagation of the linear wave, the 
computations involve a solvable Abel differential equation of the second kind. 
In the both cases we remark that the obtained solutions are not closed curves. 
This result is in the line with the recent results obtained for capillary-gravity 
water waves by using phase-plane considerations for the nonlinear system de- 
scribing the particle motion (see [TB], [II])- By the same method see also the 
results obtained for gravity water waves in [5], [10] and for constant vorticity 
gravity water waves in [T3], [TS]. Beside the phase-plane analysis, the exact so- 
lutions allow a better understanding of the dynamics (see [22], [23]). The same 
type of results are obtained for the governing equations without linearization, by 
analyzing a free boundary problem for harmonic functions in a planar domain 
(see [1] for Stokes waves, [!Ij for solitary waves and [T7] for deep-water Stokes 
waves) or by applying local bifurcation theory (see |33| for small-amplitude 
waves with vorticity). 

The existence of regular periodic travelling waves with vorticity was recently 
established (see [TT], [35] ). For steady periodic gravity waves the symmetry is 
known to be ubiquitous (see [BJ, [2T]). The study of the symmetry of rotational 
water waves was initiated in the papers [7], [E]; for irrotational flows see also 
|28j . However, exact information about the flow beneath such waves, is not 
readily available even in the irrotational case. This paper addresses this issue. 

2 Small-amplitude approximation of the water-wave problem 

The water flow under consideration is two-dimensional, bounded by a rigid hor- 
izontal surface below at z = and a free surface above at z = ho + rj(x,t), 
with ho > a constant. The undisturbed water surface is z = ho. Let 
(u(x, z,t),v(x, z,t)) be the velocity of the water and p{x,z,t) be the pressure. 
Water can be assumed to be inviscid fluid, even though it is slightly viscous. In 
problems of water waves it is also reasonable to assume that the fluid is incom- 
pressible (constant density p) (|27|). which implies the equation of mass con- 
servation (MC). A capillary-gravity wave is influenced by the effects of surface 
tension and gravity, as well as by the fluid inertia. The surface tension will play 
a role in the formulation of the boundary conditions but not in the equations of 
motion valid in the fluid domain. For the capillary-gravity water waves, the ap- 
propriate equations of motion are Euler's equations (EE)([2I]). The boundary 
conditions for the water wave problem are the kinematic boundary conditions 
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as well as the dynamic boundary condition. The kinematic boundary condi- 
tions (KBC) express the fact that the same particles always form the free water 
surface and that the fluid is assumed to be bounded below by a hard horizontal 
bed z = 0. The dynamic boundary condition (DBC) express the fact that the 
difference of pressure on the two sides of the surface 77 is balanced by the ef- 
fects of surface tension. Thus, the boundary value problem for capillary-gravity 
water waves is: 



u t + uu x + vu z = -~p x 



i" (EE) 
v t + uv x + vv z = -p> g - g 

u x +v z =0 (MC) .jv 

v = rjt + ur) x on z = h + r)(x, t) (KBC) 

v = on z = 

P = Pa - Tj , on z = ho + rj(x, t) (DBC) 

where g is the constant gravitational acceleration, po is the constant atmospheric 
pressure, the parameter T(> 0) is the coefficient of surface tension and j=r is the 
mean curvature (up to a factor 1/2) of the surface. For the surface defined as a 
function T](x, t), the mean curvature has the following expression 

R (1 + 772)3/2 U 

In respect of the well-posedness for the initial- value problem for ([1]) there has 
been significant recent progress, see [H] and the references therein. 
A key quantity in fluid dynamics is the curl of the velocity field, called vorticity. 
For two-dimensional flows we denote the scalar vorticity of the flow by 

lj(x, z)=u z — v x (3) 

In what follows we consider a flow which is uniform with depth, that is, described 
by a zero vorticity (irrotational case). 

We search for a linear approximation of the water-wave problem ([1} . First 
the system |T]) is non-dimensionalized by making use of the following scales: 
the undisturbed depth of water ho, as the vertical scale, a typical wavelength 
A, as the horizontal scale, and \fgho a s the scale of the horizontal component 
of the velocity. The surface wave itself leads to the introduction of a typical 
amplitude of the wave a. For more details see [Mj. Thus, we define the set of 
non-dimensional variables 

x 1 Y \x, z 1 Y hnz, ri 1— > aw, 1 1— > , , t, 

' ' " V¥h> ' (4) 

u 1 y sfghou, v i-> ho v/ ^ n v 

where, to avoid new notations, we have used the same symbols for the non- 
dimensional variables x, z, <q, t, u, v, on the right-hand side. 
We set the constant water density p = 1 and let us now define the non- 
dimensional pressure. If the water would be stationary, that is, u = v = 0, 
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from the equations (EE) and (DBC) with rj = 0, T = 0, we get for a non- 
dimcnsionalised z, the hydrostatic pressure po + gho(l — z). Thus, the non- 
dimensional pressure is defined by 



P ^ Po + gho(l - z) + gh p 



(5) 



Taking into account Q and ([5]) the two-dimensional capillary-gravity waves 
on water of finite depth are described, in non-dimensional variables, by the 
following boundary value problem 



u t + uu x + vu z = -p x 
8 2 (v t + uv x +vv z ) = -p z 
u x + v z = 
_ v = e(r] t + UT} X ) 
r \ n*x 

g>? J (l + £ 2«52^2)3/2 

v = 



p = e 



on z = 1 + £7?(a^, t) 
on z = 1 + £i]( x y t) 
on z = 



(6) 



where we have introduced the amplitude parameter e = and the shallowness 
parameter S = ^j-. 

For irrotational flows the vorticity equation ([3]) writes in non-dimensional vari- 
ables (0} as 

u z = S 2 v x (7) 

We observe now that, on z = 1 + erj, both v and p are proportional to e. Thus, 
with the following scaling of the non-dimensional variables, (avoiding again the 
introduction of a new notation), 



p i y ep, (u, v) i y e(u, v) 



(8) 



the problem ^ becomes 



Ut + e{uu x + vu z ) 



~Px 
-Pz 



8 2 [v t + e(uv x + vv z )] 
u x + v z = 
v = rj t + eur] x 

u = 



on z = 1 - 
on z = 1 - 
on z 



ery(a;, i) 
er)(x,t) 
= 



(9) 



and the equation ([7]) keeps the same form. Therefore, the system which describes 
the full problem in the irrotational case is given by (|H]) + ([7]). 
It is conventional to write — ^ = S 2 W e , with W e — — tt a Weber number. 
This parameter is used to measure the size of the surface tension contribution. 
By letting e — > 0, 6 and W e being fixed, we obtain a linear approximation of 
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the scaled version ©+([71) of our problem, that is, 

Ut +Px =0 
5 2 v t +p z = 
u x + v z = 

u z - 5 2 v x = (10) 
v = rjt on z = 1 

p = T) - 5 2 W e r) xx on z = 1 
v = on 2 = 

From the first three equations in (TTU)) , we get that 

v zz t = -u xz t = p xxz = -S 2 v xxt (11) 

Therefore, 

v zzt + 5 2 v xxt = (12) 

and thus 

v zz + S 2 v xx = f(x,z) (13) 

where / is an arbitrary function. Taking into account the forth equation in (llOp , 
we obtain that 

S 2 v xx = u zx (14) 
Introducing (fT4)) into ((T3|) . we have 

(v z +u x ) z = f(x, z) (15) 

and in view of the third equation in (|10p. we get that 

f(x,z)=0 (16) 

The equation (fl"3|) becomes 

v zz + S 2 v xx = (17) 

We apply the method of separation of variables, seeking a solution of this equa- 
tion in the form 

v(x,z,t) = F(x,t)G(z,t) (18) 
Substituting (TT5|) into the equation ((T7|) . we find 

rflC rfl TP 



thus, 



1 9 2 G r2 1 d 2 F 



G dz 2 5 FM (20) 
We observe in the above equation that the left hand side does not depend on z 
and the right hand side does not depend on x. Therefore, each side must be a 
constant, say 

F dx 2 ~ k ' G dz 2 ~ ^1) 
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where k > is a constant that might depend on time. With the above choice, 
the solutions of the equations in (|2"Tj) are 



F(x,t) = A sin(fcc) + B cos(fcic) 

G(x,t) = Ce kSz + De- kSz , (22) 

where A, B, C, D are constants depending on time. We made this choice of 
the sign of the constant in the equations (|21[) . in order to obtain this wave- like 
solution (|2"2")l propagating in the ir-direction. On the bed z = 0, by the last 
equation in (fTU|) . we have v = 0, thus C = —D. Therefore, 

v(x, z, t) = smh(kSz) (A sin(kx) + B coa(kx)) (23) 

where we introduced C = —D into the constants A and B. 

Taking now into the account the fifth equation in (|10j) . on z = 1 we get 

sinh(fcf)) (Asin(kx) + Bcos(kx)) = r\ t (24) 

which implies 

(^sin(fca;) + Bcos(kx)) = . J 7 * rN (25) 

smh(fco) 

Hence, 

v(x, z, t) = . sinh(fc<5z)?7 t (26) 

smh(fco) 

For the component u of the velocity field, taking into account (|26l) and the 
fourth equation of the system ([T0|) . we obtain 



u(x, z, t) = . cosh(/c5z)?7 ta; + Fix, t) (27) 

k smh(/ed) 

where F(x, t) is an arbitrary function. The components u and v of the velocity 
have to fulfill also the third equation in (fTUl) . hence, in view of (pMf and (|2"T|) , 

A , /. r x dT(x,t) kS .... , 

cosh(/cdz)?7 te2 ; H = — , cosh(kdz)rit (28) 



ksmh(kS) ' dx sinh(fc(5) 

The above relation must hold for all values of x € R, and < z < 1. It follows 

d ^ l) = (29) 
ox 

and 

Vtxx + fc 2 77t = (30) 
We seek periodic travelling wave solutions, thus, for the equation (|3"0)) with 

k = 2tt (31) 

we choose the following solution 

rjix, t) = cos(27r(a; - ct)) (32) 
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where c represents the non-dimensional speed of propagation of the linear wave 
and is to be determined. 

From (j2T))) the function T(x, t) is independent of x, therefore we will denote this 
function by Fit). 

We return now to the systems (|10p in order to find the the expressions of the 
pressure. Taking into account the first two equations in (fit)]) and the expressions 
of the velocity field from above, we obtain 

2ir5c 2 

p(x, z, t) = — — — cosh(27T(5z) cos(27r(a; — ct)) + xF'(t) (33) 

sinh(27ro) 

On the free surface z = 1 the pressure (|3"3")l has to fulfill the sixth equation of 
the system (JTDJ). Hence, in view of (|3"2"j) . we get 

2tt<5c 2 coth(27r<5) cos(27r(a; - ct)) + xT'(t) = (1 + Air 2 S 2 W e ) cos(2tt(x - ct)) (34) 

The above relation must hold for all values x € R, therefore, we get 

J-(t) = constant := cq (35) 

and we provide the non-dimensional speed of the linear wave 

2 tanh(27r<5) „ A Sr2m , A / 4yr 2 r\ , / ' 2nh \ 

c = + ^ s w ^ =^hX l+ ^) tanh br ) (36) 

We observe thus, that the speed of propagation of the wave varies with the 
wavelength A, with the undisturbed depth ho and with the coefficient of surface 
tension T. 

Summing up, the solution of the linear system (|10p is 



7](x, t) — cos(27r(a; — ct)) 

P( x > z > *) = sinhfaCa) cosh(27r,5z) cos(27r(a; - ct)) 
M ( x > z -> l ) = S inh(2^) cosh(27T(5z) cos(27r(ir - ct)) + c 
v(x, z, t) = sin l^ s) sinh(27r52;) sin(27r(a; - ct)) 



(37) 



with c given by (|36p . 



3 Exact solutions to the nonlinear equations of the motion of 
fluid particles 

Let (x(t), z(t)) be the path of a particle in the fluid domain, with location 
(x(0), z(0)) :— (xq,zq) at time t = 0. Taking into account (|37|) . the motion 
of the particle is described by the following system of nonlinear differential 
equations 

§ = u(x, z, t) = J^ S) cosh(27n5z) cos(2^(x - ct)) + c 
f = z, t) = sinh(27rfc) sm(2*(x ~ ct)) [6S) 
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The right-hand side of the differential system (|3"5j) is smooth and bounded, 
therefore, the unique solution of the Cauchy problem with initial data (xq,zq) 
is defined globally in time. 

Notice that the constant cq is the average of the horizontal fluid velocity over 
any horizontal segment of length 1, that is, 



1 



x+l 



Co = — / u(s,z,t)ds, (39) 

Jx 

representing therefore the strength of the underlying uniform current. Thus, 
Co = will correspond to a region of still water with no underlying current, 
Co > will characterize a favorable uniform current and Co < will characterize 
an adverse uniform current. 

To study the exact solution of the system (l38l) it is more convenient to re- write 
it in the following moving frame 

X = 2n(x - ct), Z = 2n5z (40) 

This transformation yields 

= sfe cosh(Z) cos(X) + 2n(c - c) 
:iI ig^sinh(Z)sin(X) 

I) c = c 

In this case, differentiating with respect to t, the system (|41[) can be written in 
the following form: 

r *x = _ sin(2X) 

\ = sinh{2 z) (42) 

^ at 2 sinh^(2ir<5) v ' 

This system integrates to 

I (f) 2 = i^S)- S h(2Z) + c 2 

ci, C2 being the integration constants. 
For the first equation in (I43I) we use the substitution 



(43) 



tan(X) = y , cos(2X) = TT ^ , sin(2X) = — ^ , dX = —^dy (44) 



— , sm(2X) = — — 
i + y l + y 1 + y 

In the new variable, the first equation in f4"5")) takes the form 



2 



dt J sinh 2 (27rr5) 



-(l- ? /)+ Cl (l + ^ (45) 



We denote by 

8ttW 



« 2 



sinh 2 (27r<5) 



(46) 
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The solution of the equation (14"51 involves an elliptic integral of first kind: 

± f- dV =t (47) 

The elliptic integral of first kind from (|47|) may by reduced to Legendre's normal 
form. In order to do this we consider first the substitution 

V 2 = s (48) 

Therefore, the left hand side in (|47|) becomes 

, <jv =± f ds 



>/(ci - a 2 )?/ 4 + 2ciy 2 + a + a 2 



(49) 

Further, we introduce a new variable ip. The dehnition of this variable depends 
on the sign of c\ — a 2 and c± + a 2 . There are three possibilities: 
ci - a 2 > 0, 

Ci — a 2 < and C\ + a 2 > 0, 
c\ — a 2 < and c\ + a 2 < 0. 

We present below only the second case, the investigation of the others will be 

presented in a future paper. 

If 

ci - a 2 < and a + a 2 > (50) 
then, we introduce the variable <p by (see [3T] Ch. VI, §4, page 602) 

a 2 + ci 2 

s = cos ip (51) 

a A — ci 

and we get 

/ 2 W , ci + a 2 \ 2a 2 (a 2 + ci) 2 . 2 2 f 2 . 2 ] 

(ci - a )s(s + 1) s H ^ = — r-5 sin 93 cos w 1 — fc, sin y> 

\ ci — a z / (a z — ci)^ 

2(a 2 + Cl ) . 

as = sin (p cos yayj 

a 2 — c\ 

where the constant < k\ < 1 is given by 

Therefore we obtain the Legendre normal form of the integral in (|47p. that is, 

1 

V2a 



±4r~ l , ^ (53) 



1 — k{ sin ip 
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The inverse of the integral in (|53p is the Jacobian elliptic function sine amplitude 
(see, for example, pQ), an odd periodic function of order two, 



sn [±V2at; kij := sinip (54) 
In view of the notations (|48p. (jSTj) . we get that 



y(t) = ±y ^ cn (± V2a i; fc x ) = ±y ^ cn (\/2a t; fci) (55) 

cn being the Jacobian elliptic function cosine amplitude, an even periodic func- 
tion of order two. 

For the second equation in (|4"5|) we use the substitution 

tanh(Z) = w , cosh(2Z) = 1 + , dX = — (56) 

1 — w 2 1 — ur 

In the new variable, the second equation in (1431) takes the form 



dw^ 2 



-(l-w A ) + c 2 (l-w 2 ) 2 (57) 



dt J sinh 2 (27r<S) 
The solution of the equation (|57|) involves an elliptic integral of first kind: 



±[- dW =t (58) 

J V ( c 2 - a 2 )w 4 - 2c 2 w 2 + c 2 + a 2 

where a 2 is the constant from (|46|) . The elliptic integral of first kind from (|58[) 
may by reduced to Legendre's normal form. In order to do this we consider first 
the substitution 



w 2 = r (59) 



The left hand side in (l58l) becomes 



dw f ds 



2^/{c 2 - a 2 )w 4 - 2c 2 u; 2 + c 2 + a 2 J 2 ^(c 2 ~ a 2 )r(r - l)(r - g±^) 

C2 (60) 

As in the case of the integral in (|49|) . we introduce a new variable <fi. The 
definition of (f> depends on the sign of c 2 — a 2 and c 2 + a 2 . There are three 
possibilities: 

c 2 - a 2 > 0, 

c 2 - a 2 < and c 2 + a 2 > 0, 
c 2 — a 2 < and c 2 + a 2 < 0. 

We present below only the last case, the investigation of the others will be 

presented in detail in a future paper. 

If 

c 2 - a 2 < and c 2 + a 2 < (61) 
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then, we introduce the variable <f> by (see [3T] Ch. VI, §4, page 602) 

2a 2 . o , 



r = 1 



c 2 - a z 



■ sin 



(62) 



and we get 



(c 2 - a )r{r - 1) I r - 



c 2 + a' 
c 2 - a 2 



4a 4 



(a 2 - c 2 ) 



sin 2 6 cos 2 i 



(l - k\ sin 2 



4a 2 



c 2 - a" 

where the constant < fc 2 < 1 is given by 

2a 2 



k l - 



(63) 



a — c 2 

Therefore we obtain the Legendre normal form of the integral in ([58)1. that is, 

1 



\/a 2 - c 2 

The inverse of the integral in (|64[) is 



1 - k\ sin 2 



f 



sn (± y/ a 2 - c 2 t; k 2 
In view of the notations ([55)) . (|6"2")l . we get that 



t«(t) =±4/1 



2a 2 



a 2 — c 2 



(64) 



(65) 



(66) 



Therefore, from (|4"4")) and (|56l) . the solution of the system (1431) has the following 
expression 

= arctan [#(*)] 



Z(t) = arctanh [tu(t)] = \ In ±±g^ 



(67) 



(*) 



with y{t) given by (|55l) and w(t) given by ([66| . From (|40)) and (l67l) . the solution 
of the system (j3HJ) with the constant Co equals the speed of propagation of the 
linear wave c, have the following expressions: 



x(t) = ct± ±arctan 
z(t) = ±2^5 arctanh 



££±£Lcn 

a^ — ci 



(V2at;fci) 



1-^ sn 2 (V^^;fc 2 ) 



(68) 



We remark that the curve in ([68)) is not a closed curve. This result is in the line 
with the results obtained in g], [5J, [ID], [H], [15], [17], [18], [19], [22], [23]. 
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II) c ^ c 

Differentiating with respect to t the system (|4Tj) we get 



d 2 X dX 

-p- +6tanp0— + a 2 sin(2X) - 6 2 tan(X) = (69) 

where a 2 is the constant from (j46l ) and 

b := 27r(c - c) (70) 
Using the substitution (|44[) . the equation (|69|l takes the form 

<r-u 2y fdy\\ by dji + 2a2y _ bMl + y2) ^ (71) 



di 2 l + y 2 \dt J dt 
For (see [5S], 6.54, page 554) 

*»-$, (72) 
the equation ([7T|) becomes an Abel differential equation of the second kind 

= txV - 5yp - 2a2y + ^ + y 2 ^ w 

ay 1 + y 2 
The substitution (see [25]. 4.11. pages 26-27) 

u (y) = p(y)E(y), where E(y) = exp f- J - 2y 

brings this equation to the simpler form 



d V = i I 2 ' ( 74 ) 



y / i + v 



u ^ = - b ^— u -2a 2 - y ^ 2 + b 2 -^ (75) 
dy l + y 2 (l + y 2 ) 2 i + y 2 



The equation (|75|) with the substitution 



i + y 2 7 y 2 

can be written in the canonical form: 

d« 2a 2 /2A 

exp UJ (77) 

The equation (|77|) is solvable (see [29 , 8., page 111), its solution can be written 
out in the following parametric form 
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C being a constant. From (fTB"]) and ([75)1 we get the expression of y 

y(r)=±< 



t 2 - 2a? 



(C - Mn|r + Vr 2 - 2a 2 \)' 



(79) 



By (|72]) and (US), we have J£ = (f;) and therefore, from dZSJ) and f79]). 

the relation between t and r is the following: 



< = 



Vt 2 - 2a? ^t 2 - 2a 2 ~ (C - bin |r + Vt 2 - 2a 2 |) 2 
Thus, taking into account (|4"4"]) . we obtain 

X(t) = arctan [y(t)], 



(80) 



(81) 



with y(r) given by (1791 and r given implicitly by (|5U|) . 

In order to determine Z(i) from the system (|4Tj) . with (f8T|) in view, we write 
the second equation of this system in the form 



dZ 



sinh(Z) sinh(27r5) 
Integrating, we get 



4tt 2 <5c . , 4ir 2 5c y(t) 

sm(arctan [y(t)\) at = — 



In 



tanh 



3inh(27T(5) y/l + y*(t) 
An 2 dc y(t) 



If 



then 



sinh(27T(5) ^/l + y 2 (t) 
4tt 2 5c y(t) 



dt + const 



dt (82) 



(83) 



Z(t) = 2arctanh 



sinh(27T<5) ^/l + y 2 ^) 
47T 2 (5c 



(it + const < 



(84) 



exp 



y(t) 



■. dt + const 



sinh(27r5) ^l + y 2 {t) 
From P0|) , (|81l) and (1551) . the solution of the system ([35]) is written now as 
= ct + 7^: arctan \y(t)} 



(85) 



z(i) = -^arctanh 



exp ( f . ^ (t) , + const ) 



(86) 



with y{r) given by (f79|) and r given implicitly by (|80|). 

We remark that the curve in ([86]) is not a closed curve. This result is in the line 
with the results obtained in g], [5], [10], Q3], [23, PU, US, PS], [21], [23]. 
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